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Abstract 

We study bifurcations of periodic travelling waves in granular dimer chains from the anti- 
continuum limit, when the mass ratio between the light and heavy beads is zero. We show that 
every limiting periodic wave is uniquely continued with respect to the mass ratio parameter 
and the periodic waves with the wavelength larger than a certain critical value are spectrally 
stable. Numerical computations are developed to study how this solution family is continued 
to the limit of equal mass ratio between the beads, where periodic travelling waves of granular 
monomer chains exist. 
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1 Introduction 

Wave propagation in granular crystals has been studied quite intensively in the past ten years. 
Granular crystals are thought to be closely-packed chains of elastically interacting particles, which 
obey the Fermi-Pasta-Ulam (FPU) lattice equations with Hertzian interaction forces. Experimen- 
tal work with granular crystals and their numerous applications [6l [21] stimulated theoretical and 
mathematical research on the granular chains of particles. 

Existence of solitary waves in granular chains was considered with a number of analytical and 
numerical techniques. In his two-page note, MacKay [IB] showed how to adopt the technique of 
Priesecke and Wattis [8j to the proof of existence of solitary waves. English and Pego [7] used 
these results to prove the double-exponential decay of spatial tails of solitary waves. Numerical 
convergence to the solitary wave solutions was studied by Ahnert and Pikovsky [1]. Stefanov and 
Kevrekidis [23] reviewed the variational technique of [8] and proved that the solitary waves are 
bell-shaped (single-humped) . 

Recently, the interest to granular crystals has shifted towards periodic travelling waves as well 
as travelling waves in heterogeneous (dimer) chains, as more relevant for physical experiments 
[U [20]. Periodic wave solutions of the differential advance-delay equation were considered by 
James in the context of Newton's cradle [11] and homogeneous granular crystals [12]. In par- 
ticular, numerical approximations in [12] suggested that periodic waves with wavelength larger 
than a critical value are spectrally unstable. Convergence to solitary waves in the limit of infinite 
wavelengths was also illustrated numerically and asymptotically in [12]. More recent work |13j 
showed non-existence of time-periodic breathers in homogeneous granular crystals and existence 
of these breathers in Newton's cradle, where a discrete p-Schrodinger equation provides a robust 
approximation. 

Periodic waves in a chain of finitely many beads closed in a periodic loop were approximated 
by Starosvetsky et al. in monomers [22] and dimers |14j by using numerical techniques based 
on Poincare maps. Interesting enough, solitary waves were found in the limit of zero mass ratio 
between lighter and heavy beads in [14]. It is explained in [14] that these solitary waves are 
in resonance with linear waves and hence they do not persist with respect to the mass ratio 
parameter. Numerical results of [H] indicate the existence of a countable set of the mass ratio 
parameter values, for which solitary waves should exist, but no rigorous studies of this problem 
have been developed so far. Recent work [15] contains numerical results on existence of periodic 
travelling waves in granular dimer chains. 

In our present work, we rely on the anti-continuum limit of the FPU lattice, which was recently 
explored in the context of existence and stability of discrete multi-site breathers by Yoshimura 
|24| . By using a variant of the Implicit Function Theorem, we prove that every limiting periodic 
wave is uniquely continued with respect to the mass ratio parameter. By the perturbation theory 
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arguments (which are similar to the recent work in |19j). we also show that the periodic waves with 
the wavelength larger than a certain critical value are spectrally stable. Our results are different 
from the asymptotic calculations in [T3|, where a different limiting solution is considered in the 
anti-continuum limit. 

The family of periodic nonlinear waves bifurcating from the anti-continuum limit are shown 
numerically to extend all way to the limit of equal masses for the dimer beads. The periodic 
travelling waves of the homogeneous (monomer) chains considered in jl2j are different from the 
periodic waves extended here from the anti-continuum limit. In other words, the periodic waves 
in dimers do not satisfy the reductions to the periodic waves in monomers even if the mass ratio is 
one. Similar travelling waves consisting of binary oscillations in monomer chains were considered 
a while ago with the center manifold reduction methods |10j . 

The paper is organized as follows. Section 2 introduces the model and sets up the scene for the 
search of periodic travelling waves. Continuation from the anti-continuum limit is developed in 
Section 3. Section 4 gives perturbative results that characterize Floquet multipliers in the spectral 
stability problem associated with the periodic waves near the anti-continuum limit. Numerical 
results are collected together in Section 5. Section 6 concludes the paper. 



2 Formalism 
2.1 The model 

We consider an infinite granular chain of spherical beads of alternating masses (a so-called dimer), 
which obey Newton's equations of motion, 

mXn = V'{yn - Xn) - V'{Xn - Vn-l), rj, /-, x 

My„ = y'(x„+i - Vn) - V'iVn - Xn), ^ ' ^ > 

where m and M are masses of light and heavy beads with coordinates {xn}nez and {yn}nez, 
respectively, whereas V is the interaction potential. The potential V represents the Hertzian 
contact forces for perfect spheres and is given by 

V{x) = ^|x|i+-F(-x), (2) 

where a = | and H is the Heaviside step function with H{x) = 1 for x > and H{x) = for 
X < 0. The mass ratio is modeled by the parameter := Using the substitution, 

reeZ: x„(t) = n2„_i(r), yn{t) = ew2n{T), t = ^/mT, (3) 

we rewrite the system of Newton's equations ([1]) in the equivalent form: 

U2n-1 = y'{£W2n - U2n-l) -V'{u2n-1 - £W2n-2), ^ ^ U) 
W2n = £y'{u2n+l- £W2n) - £V'{£W2n- U2n-l), 

The value e = correspond to the anti-continuum limit, when the heavy particles do not move. 
At the limit of equal mass ratio e = 1, we note the reduction, 

nGZ: U2n~l{T) = U2n~l{T), W2n{T) = U2n{T) , (5) 
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for which the system of two granular chains reduces to the scalar granular chain (a so-called 
monomer), 

tin = V'{Un+l - Un) " V {Un " C/n-l), n G Z. (6) 

The system of dimer equations ([4]) has two symmetries. One symmetry is the translational 
invariance of solutions with respect to r, that is, if {u2n-i{'T)-,W2n{T)}n& is a solution of (jU, then 

{u2n-l{T + h),W2n{T + h)}n& (7) 

is also a solution of ^ for any 6 € M. The other symmetry is a uniform shift of coordinates 
{u2n-i-,W2n}n(i7. in the direction of (e, 1), that is, if {n2„_i(T), ^D2n(T")}ngz is a solution of (jl]), 
then 

{u2n-l{T) + ae, W2n{T) + a}nez (8) 

is also a solution of dH for any a G M. 

The system of dimer equations (j4]) has the symplectic structure 

dH dH dH dH ^ , ^ 

U2n-1 = 7^ , P2n-1 = , W2n = 7^ , q2n = "7^ , n G (9) 

OP2n-l OU2n-l Oq2n OW2n 

where the Hamiltonian function is 

H = 1^^ + lln) + X] ^(^"^2n - U2n-l) + ^ V {u2n-l - £W2n-2), (10) 

neZ neZ neZ 

written in canonical variables {u2n-i,P2n-i = U2n-i,W2n,q2n = 'W2n}n&z- 
2.2 Periodic traveling waves 

We shall consider 27r-periodic solutions of the dimer system @, 

U2n-l{T) = U2n-l{T + 2tx), l(;2n (t) = W2n (t + 27r) , T S M, n G Z. (11) 

Travelling waves correspond to the special solution to the dimer system (jj]), which satisfies the 
following reduction, 

U2n+l{T) = U2n-l{T + 2q), W2n+2 (t) = t(;2n (t + 2g) , T G M, n G Z, (12) 

where q G [0, vr] is a free parameter. We note that the constraints (jlip and ()12p imply that there 
exists 27r-periodic functions and w.^ such that 

'U2n-i(T) = u*(t + 2(/n), W2n(''") = ''^♦(t + 2(7n), r G M, n G Z. (13) 

In this context, q is inverse proportional to the wavelength of the periodic wave over the chain 
n £ 7j. The functions and satisfy the following system of differential advance-delay equations: 

u*{t) = V'{ew^{T) - u^{t)) -V'{u^{t) - ew^ir -2q)), 

ii,^^r) = eV'{u,{T + 2q)-ew,{T))-eV'{ew,{T)-u,{T)), ^ ' 
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Remark 1. A more general traveling periodic wave can be sought in the form. 

U2n-i{T) = u^{cT + 2qn), W2n{T) = w^{cT + 2qn), T G M, n G Z, 

where c > is an arbitrary parameter. However, the parameter c can be normalized to one thanks 
to invariance of the system of dimer equations ^ with respect to a scaling transformation. 

Remark 2. For particular values q = where m and N are positive integers such that 1 < 
m< N , periodic travelling waves satisfy a system of2mN second-order differential equations that 
follows from the system of lattice differential equations ^ subject to the periodic conditions: 

U_i = U2mN~l, U2mN+l = m, Wq = W2mN , W2mN+2 = W2- (15) 

This reduction is useful for analysis of stability of periodic travelling waves and for numerical 
approximations. 



2.3 Anti-continuum limit 

Let be a solution of the nonlinear oscillator equation, 

(p = V'{-^)-V'{^) ^ (^+|(^r-V = 0. (16) 

Because a = |, bootstrapping arguments show that if there exists a classical 27r-periodic solution 
of the differential equation (fT6|) . then ip G Cpgj,(0, 27r). 

The nonlinear oscillator equation (jl6p has the first integral, 



E = \^' + -r^Wr'. (17) 

The phase portrait of the nonlinear oscillator ()16p on the (93, (/9)-plane consists of a family of closed 
orbits around the only equilibrium point (0, 0). Each orbit corresponds to the T-periodic solution 
for if, where T is determined uniquely by energy E. It is well-known \12\ [24] that, for a > 1, the 
period T is a monotonically decreasing function of E such that T — t- 00 as £^ — t- and T — t- as 
E ^ 00. Therefore, there exists a unique Eq G M+ such that T = 2tt for this E = Eq. We also 
know that the nonlinear oscillator (jl6p is non-degenerate in the sense that T'{Eq) ^ (to be more 
precise, T'{Eo) < 0). 

In what follows, we only consider 27r-periodic functions (p which are defined by (I17p for E = Eq. 
For uniqueness arguments, we shall consider initial conditions (p{0) = and ^p^O) > 0, which 
determine uniquely one of the two odd 27r-periodic functions ip. 

The limiting 27r-periodic travelling wave solution at e = should satisfy the constraints (|12p . 
which we do by choosing for any fixed q G [0, vr], 

e = 0: U2n~i{T) = ip{T + 2qn), w2n{T) = 0, tGM, nGZ. (18) 

To prove the persistence of this limiting solution in powers of e within the granular dimer chain 
dH), we shall work in the Sobolev spaces of odd 27r-periodic functions for {u2n~i}n£Z, 

H^ = ^u£ -f/per(0, 27r) : u{-t) = -u{t), r G m} , G No, (19) 
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and in the Sobolev spaces of 27r-periodic functions with zero mean for {w2n}nez, 

H'^ = l^we i/pcr(0, 27r) : w{T)dT = o| , A; G Nq. (20) 

The constraints in (I19p and ()20p reflects the presence of two symmetries ([7]) and ([8]). The two 
symmetries generate a two-dimensional kernel of the linearized operators. Under the constraints 
in (jl9p and (j20p . the kernel of the linearized operators is trivial, zero-dimensional. 

It will be clear from analysis that the vector space defined by ([20]) is not precise enough to 
prove the persistence of travelling wave solutions satisfying the constraints (|12p . Instead of this 
space, for any fixed q G [0, vr], we introduce the vector space by 

Hi = [weHl^^{Q,2T:): ^r) = -u;(-r - 2g)} , A: G Nq. (21) 

We note that C H^, because if the constraint w{t) = —w{—t — 2q) is satisfied, then the 
27r-periodic function w has zero mean. 



2.4 Special periodic traveling waves 

Before developing persistence analysis, we shall point out three remarkable explicit periodic trav- 
elling solutions of the granular dimer chain 1^ for q = 0, q = ^ and q = ir. For q = |, we have an 
exact solution 

vr 

^=2' '^2n^l{T) = ip{T + 7^71), W2n{T) = 0. (22) 

This solution preserves the constraint V'{u2n+i) = ^'(— 'U2n-i) in equations (jH thanks to the 
symmetry (p{T — vr) = (p{T -|- vr) = — y?(r) on the 27r-periodic solution of the nonlinear oscillator 
equation (fT6|) . 

For either q' = or q' = vr, we obtain another exact solution, 

g = {0,7r}: U2n-i{r) = jr^^^^^, W2n(T) = ^^^^2^3 ' (23) 

By construction, these solutions (j22p and (j23p persist for any e > 0. We shall investigate 
if the continuations are unique near e = for these special values of q and if there is a unique 
continuation of the general limiting solution (jlSp in e for any other fixed value of g G [0, vr]. 

Furthermore, we note that the exact solution (I23p for q = vr at e = 1 satisfies the constraint ([5]) 
with U2n-i{T) = —U2n{T) = U2n{'T — T^)- This reduction indicates that the function (j23p for e = 1 
satisfies the granular monomer chain ([6]) and coincides with the solution considered by James |12j . 
On the other hand, the exact solutions (j22p for q = f and (j23p for g = do not produce any 
solutions of the monomer chain at e = 1. This indicates that there exists generally two distinct 
solutions at e = 1, one is continued from e = and the other one is constructed from the solution 
of the monomer chain ([6]) in [12] . 
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3 Persistence of periodic traveling waves near e = 
3.1 Main result 

We consider the system of differential advance-delay equations (jl4p . The limiting solution (jlSp 
becomes now 

e = 0: u^{t) = ^{t), w^{t) = 0, r G M, (24) 

where (pis a unique odd 27r-periodic solution of the nonlinear oscillator equation (fT6]) with 0(0) > 0. 
We now formulate the main result of this section. 

Theorem 1. Fix q € [0, vr]. There is a unique continuation of 2tt -periodic traveling wave \24^ 
in £, that is, there is a Sq > such that for every e G (0, eo); there are C > and a unique 
2iT-periodic solution {u^^w^,) G x H"^ of the system of differential advance-delay equations (T^ 
such that 

ll^i* — fpll ff2 < Ce^, ||w;*||t72 < Ce. (25) 



Remark 3. By Theorem\^ the limiting solution ^24\ ) for q G |0, §,7r} is uniquely continued in 
£. These continuations coincide with the exact solutions I122\) and 

3.2 Formal expansions in powers of e 

Let us first consider formal expansions in powers of e to understand the persistence analysis from 
£ = 0. Expanding the solution of the system of differential advance-delay equations (|14p . we write 



u,(r) = 99(T)+e2-ul^^(r) + o(e2), w,{t) = ewi^\T) + o{e^), (26) 
and obtain the linear inhomogeneous equations 

wI}\t) = F(i)(t) := V\^{t + 2q)) - V'{-^{t)) (27) 

and 

#(t) + a\^{TT-^uf\T) = Fi^\T) := V" {-ip{T))wi'\T) + V" {ip{T))wi'\T - 2q). (28) 

Because V is but not , we have to truncate the formal expansion (j26p at o(e^) to indicate 
that there are obstacles to continue the power series beyond terms of the 0{£^) order. 
Let us consider two differential operators 

Lo = -^ ■■ //^^,(0,27r)^L2^,(0,27r), (29) 

L = + a\^{T)r-' : H^^M^tt) ^ lI^M^tt), (30) 

As a consequence of two symmetries, these operators are not invertible because they admit one- 
dimensional kernels, 

Ker(Lo) = span{l}, Ker(L) = span{(/3}. (31) 
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Note that the kernel of L is one-dimensional under the constraint T'(Eq) ^ (see Lemma 3 in 

|12j for a review of this standard result). 

To find uniquely solutions of the inhomogeneous equations (|27p and (j28p in function spaces 
and respectively, see (fT9]) and (f20l) for definition of function spaces, the source terms must 

satisfy the Fredholm conditions 

(l,FW)i2^^=0 and {ip,Fi^'))L2^^=0. 
The first Fredholm condition is satisfied, 

r>27r f2n r^ir 



/ [V'i^ir + 2g)) - V'i-^ir))] dr = / Vi^pir + 2g))dr - / V 
Jo Jo Jo 



{-^{T))dT = 0, 



because the mean value of a periodic function is independent on the limits of integration and 
the function (/9 is odd in r. Since F^^ e L2 , there is a unique solution w^^^ G H"^ of the linear 
inhomogeneous equation ()27p . 

The second Fredholm condition is satisfied, 



0(r) ry"(-(^(r))4'^ (r) + y"(v9(r))4'^ (r - 2g)' 







dr = 0, 



if the function Fu'^ is odd in r. If this is the case, then Fu'^ S L\ and there is a unique solution 
^(2) g ^2 q£ linear inhomogeneous equation (j28p . To show that -Pi^^ is odd in r, we will prove 
that w')^'^ satisfies the reduction 

4i)(r) = -t/;«(-r-2g), ^ )(-r) = -Fi2)(r), r € M. (32) 

It follows from the linear inhomogeneous equation ()27p that 

ti)i'^(T) + w^}\-T - 2q) = V'{^{t + 2q)) - V'{-^{t)) + V'{^{-t)) - V'{-ip{-T - 2q)) = 0, 

where the last equality appears because is odd in r. Integrating this equation twice and using 
the fact that w^^^ G H^, we obtain reduction ()32p . Note that the reduction (I32p implies that 
wi^^ G ^2 , where C i^^ given by (j2T]) . 

3.3 Proof of Theorem [1] 

To prove Theorem [U we shall consider the vector fields of the system of differential advance-delay 
equations (fH|) . 

FMr),w{T),e) := V'iewir) - u{r)) - V'{u{t) - ew{T - 2q)), 

FUu{T),w{T),e) := eV'iu{T + 2q) - ew{T)) - eV'{ew{T) - n(r)), ^ ^ ^"'''^ 
We are looking for a strong solution {u^:,w^:) of the system (jl4p satisfying the reduction, 

n*(— r) = — n*(T), w^,{t) = -w^,{—t — 2q), r G M, (34) 
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that is, € Hl{R) and £ Hl{R). 

If {u,w) € //^ X iJ^, then is odd in r. Furthermore, since y is C^, then is a map 
from if^ X x M to and its Jacobian at e = is given by 

DuFu{u,w,0) = V"{-u) -V"{u) = -a\u\°'-^, u-, 0) = 0. (35) 

On the other hand, under the constraints (I34p . we have F^ G L^, because 



/ F^(n(r),u;(T),e)(ir = £ V'{u{t + 2q) + £w{-t - 2q))dT 
Jo Jo 



2n 

V'{£w{t) + u{-T))dT = 0. 

Moreover, under the constraints (j3l|) . we actuaUy have F^ S because 

F^{u{t),w{t),£) + F^(n(-T - 2g),?i;(-r - 2g),e) 

= ey' (ti(T + 2q) - £w{t)) - £V'{£w{t) - n(r)) 

+ £V'{u{-t) - £w{-T - 2q)) - £V'{£w{-T - 2q) - u{-T - 2q)) 

= 0. 

Since V Is C"^, then F^) is a map from H'^^ x to and its Jacobian at e = is given by 

DuF^{u, w, 0) = 0, D.^F^{u, w, 0) = 0. (36) 
Let us now define the nonlinear operator 



fuiu, w, e) := ^ - Fu{u, w, e) 
Uiu,w,£) := ^ - F^{u,w,£ 



We have (/„, /^) : i^^ x ^2 ^ M ^ x are near the point (99, 0, 0) £ x x M. To 
apply the Implicit Function Theorem near this point, we need {fu,fw) = at {u,w,£) = ((^,0,0) 
and the invertibility of the Jacobian operator {fu, fw) with respect to (n, w) near (</?, 0,0). 

It follows from ([35]) and ([36]) that the Jacobian operator of (/„, fw) at (99, 0, 0) is given by the 
diagonal matrix of operators L and Lq defined by (j29p and (j30p . The kernels of these operators 
in ([3T]) are zero-dimensional in the constrained vector spaces (fT9|) and ([20]) (we actually use space 
([2T]) in place of space (1201)). 

By the Implicit Function Theorem, there exists a continuation of the limiting solution (|24p 
with respect to e as the 27r-periodic solutions {u^:,w^) £ x of the system of differential 
advance-delay equations ([H]) near £ = 0. From the explicit expression ([33]) . we can see that 
11'"^* 11^2 = 0(e) whereas — (/j||/i-2 = O(e^) as e — )• 0. The proof of Theorem [1] is complete. 
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4 Spectral stability of periodic traveling waves near e = 
4.1 Linearization at the periodic traveling waves 

We shall consider the dimer chain equations ([4]), which admit for small e > the periodic trav- 
eling waves in the form (|13p . where {u^,w^,) is defined by Theorem [TJ Linearizing the system of 
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nonlinear equations (jH) at the periodic traveling waves (jl3p . we obtain the system of linearized 
dimer equations for small perturbations, 

M2„-i = V"{eW:,{T + 2qn) - u^.{t + 2qn)){ew2n - U2n-i) 

- V"{u^{t + 2qn) - sw^{t + 2qn - 2q)){u2n-i - ew2n-2), /^gx 
W2n = £y''{U:^{T + 2qn + 2q) - ew^{T + 2qn)){u2n+i - £W2n) 

- eV"{ew^.{T + 2qn) - u^.{t + 2qn)){ew2n - U2n-i)-, 

where n S Z. A technical complication is that V" is continuous but not continuous differentiable. 
This will complicate our analysis of perturbation expansions for small e > 0. Note that the 
technical complications does not occur for exact solutions (I22p and ()23p . Indeed, for exact solution 
(j22p with 9 = §, the linearized system ()38p is rewritten explicitly as 



U2n-1 + alv^r ^U2n-l = £ {y"{-^)w2n + V" {Lp)w2n-2) , /ggx 
id2n + 2e^V"{-ip)w2n = eV"{-ip){u2n+l + W2n-l). ^ ^ 



For exact solution (j23p with q = or q = n, the linearized system (j38p is rewritten explicitly as 

\ *2n + l^|</'r-W = i:^(V^"M^X2n+l+V^"(-¥^)n2n-l). ^ ^ 

In both cases, the linearized systems ()39p and (j40p are analytic in e near e = 0. 

The system of linearized equations ()38p has the same symplectic structure ([9]) , but the Hamil- 
tonian is now given by 

^ = l^{pln-i + Q2n) + l^y"{£w^i-r + ^qn)-u4T + 2qn)){ew2n-U2n-i? 

+ \Y1 ^"(M^ + 2gra) - ew^T + 2qn - 2q)){u2n-i - ew2n-2? ■ (41) 
The Hamiltonian H is quadratic in canonical variables {u2n-i,P2n-i = U2n-i,W2n,Q2n = i^2n}„gz- 



4.2 Main result 

Because coefficients of the linearized dimer equations (j38p are 27r-periodic in r, we shall look for 
an infinite-dimensional analogue of the Floquet theorem that states that all solutions of the linear 
system with 27r-periodic coefficients satisfies the reduction 

u(r + 27r) = A^u(r), r G M, (42) 

where u :=[••• , W2n-2,U2n-i,W2n, U2n+i, ' ' ' ] and is the monodromy operator. 

Remark 4. Let q = ^ for some positive integers m and N such that 1 < m < N . In this case, 
the system of dimer equations ^ can be closed into a chain of 2mN second-order differential 
equations subject to the periodic boundary conditions /[T5\). Similarly, the linearized system i T^) 
can also be closed as a system of 2mN second-order equations and the monodromy operator A4 
becomes an infinite diagonal composition of a AmN-by-AmN Floquet matrix, each matrix has AmN 
eigenvalues called the Floquet multipliers. 
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We can find eigenvalues of the monodromy operator M by looking for the set of eigenvectors 
in the form, 

U2n-l{r) = U2n~l{r)e^\ U2n(.r) = W2n(.T)e^\ TGM, (43) 

where {U2n-i,W2n) are 27r-periodic functions and the admissible values of A are found from the 
existence of these 27r-periodic functions. The admissible values of A are called the characteristic 
exponents and they define the Floquet multipliers fi by the standard formula /U = e^'^'^. 
Eigenvectors ()43p are defined as 27r-periodic solutions of the linear eigenvalue problem, 

' U2n-i + 2\U2n-i + >?U2n~i = V" {ew^{T + 2qn) - n,(T + 2qn)){eW2n - f/2n-l) 
-V"{u^{T + 2qn)-ew^{T + 2qn-2q)){U2n-i-eW2n~-2), 
W2n + 2\W2n + \^W2n = £y"{u,{T + 2qn + 2q)-ew,{T + 2qn)){U2n+i-eW2n) ^ ' 
- eV"{ew^{T + 2qn) - u^\t + 2qn)){eW2n - U2n-i)- 

The Krein signature, which plays an important role in the studies of spectral stability of 
periodic solutions (see Section4 in [2]), is defined as the sign of the 2-form associated with the 
symplectic structure ([9]): 

(J = i ^ [U2n-lp2n-l " ^271-1^271-1 + W2nq2n " W'2nQ2n] , (45) 

where {u2n-i-,P2n-i = U2n-i,'W2n,Q2n = 'W2n}nez eigenvector (1431) associated with an eigen- 
value A G i]R+. Note that by the symmetry of the linear eigenvalue problem (j44p . it follows that 
if A is an eigenvalue, then A is also an eigenvalue, whereas the 2-form a is constant with respect 
to T G M. 

If e = 0, the monodromy operator in (j42p is block-diagonal and consists of an infinite 
set of 2-by-2 Jordan blocks, because the dimer system (jH) is decoupled into a countable set of 
uncoupled second-order differential equations. As a result, the linear eigenvalue problem ()44p 
with the limiting solution (llSp admits an infinite set of 27r-periodic solutions for A = 0, 

e = 0: =C2n-ivi(T + 2gn), W^,^ = a2n, nGZ, (46) 

where {c2n-i,CL2n}nez are arbitrary coefficients. Besides eigenvectors (Ii6]l . there exists another 
countable set of generalized eigenvectors for each of the uncoupled second-order differential equa- 
tions, which contribute to 2-by-2 Jordan blocks. Each block corresponds to the double Floquet 
multiplier /i = 1 or the double characteristic exponent A = 0. When e ^ but e ^ 1, the charac- 
teristic exponent A = of a high algebraic multiplicity splits. We shall study the splitting of the 
characteristic exponents A by the perturbation arguments. 
We now formulate the main result of this section. 

Theorem 2. Fix q = ^ for some positive integers m and N such that 1 < m < N . Let 
{u^:,W:^) £ X be defined by Theorem{^for sufficiently small positive e. Consider the linear 
eigenvalue problem \^J^ subject to 2mN -periodic boundary conditions U5\} . There is a Eq > 
such that, for every e G (0,eo), there exists qo{£) G (O, ^) such that for every q G {O,qo{e)) or 
g G (vr — qo{e),7r], no values of X with Re(A) / exist, whereas for every q G (g'o(e))7r — loi^)), 
there exist some values of A with Re(A) > 0. 
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Remark 5. By Theorem\^ periodic traveling waves are spectrally stable for q S (0, go(^)) md 
q £ (tt — go(e)>7r] and unstable for q € {qo{£),7r — qoi^))- Therefore, the linearized system for 
the exact solution i22\) with 9 = f subject to ^-periodic boundary conditions (m = 1 and N = 2) 
is unstable for small e > 0, where the linearized system J^Op for the exact solution l[23\) with q = ir 
subject to 2-periodic boundary conditions (m = 1 and N = 1) is stable for small e > 0. 

Remark 6. The result of Theorem\^is expected to hold for all values of q in [0, vr] but the spectrum 
of the linear eigenvalue problem for the characteristic exponent A becomes continuous and 
connected to zero. An infinite- dimensional analogue of the perturbation theory is required to study 
eigenvalues of the monodromy operator Ai in this case. 

Remark 7. The case q = is degenerate for an application of the perturbation theory. Neverthe- 
less, we show numerically that the linearized system for the exact solution 123\} with q = 
(m = 1 and N ^ oo) is stable for small e > and all characteristic exponents are at least double 
for any e > 0. 



4.3 Formal perturbation expansions 

We would normally expect splitting A = if the limiting linear eigenvalue problem at 

e = is diagonally decomposed into 2-by-2 Jordan blocks [19]. However, in the linearized dimer 
problem this splitting occurs in a higher order, that is, A = 0{e), because the coupling 

between the particles of equal masses shows up at the 0{e^) order of the perturbation theory. 
Regular perturbation computations in O(e^) would require V" to be at least C^, which we do 
not have. In the computations below, we neglect this discrepancy, which is valid at least for 
q = and q = tt. For other values of q, the formal perturbation expansion is justified with the 
renormalization technique (Section 4.6). 

We expand 27r-periodic solutions of the linear eigenvalue problem (|44p into power series of e: 

A = eA«+e2;^(2)+o(e2) (47) 

and 

C/2n-i = <li + eUiZi + Mil + o{e% .... 

where the zeroth-order terms are given by (j46|) . To determine corrections of the power series 
expansions uniquely, we shall require that 

{^Mn-i)Ll^, = {hwg^)L2^^ = 0, neZ, i = l,2. (49) 

(7) 

Indeed, if C^2n-i contains a component, which is parallel to if, then the corresponding term only 
changes the value of C2n-i in the eigenvector (|45|) . which is yet to be determined. Similarly, if a 

(7) (7) 
27r-periodic function IV2n ^ nonzero mean value, then the mean value of only changes 

the value of a2n in the eigenvector ([I6l) . which is yet to be determined. 
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The linear equations (|44p are satisfied at the 0{e^) order. Collecting terms at the 0{e) order, 
we obtain 

+ V"{-^{t + 2qn))wi^ + V"{^{t + 2qn))wi^_„ (50) 
[ = -2A(i)#£) + V"{^{t + 2qn + 2g))C/£Vi + V"{-^{t + 2^n))C/(li. 

Let us define solutions of the following linear inhomogeneous equations: 

v + a\ip\'^-'^v = -2^, (51) 
y± + a\^r~'y± = V"{±^), (52) 
z± = V"{±Lp)ip. (53) 

If we can find uniquely 27r-periodic solutions of these equations such that 
then the perturbation equations ([50]) at the 0{e) order are satisfied with 



^2n-i = C2n-iA(^)f (r + 2qn) + a2„y-.(r + 2qn) + a2n~2y+{r + 2q'n), 
^2n = C2n+i^+(T + 2gn + 2q) + C2„_iz_(t + 2qn). 



(54) 



The linear equations (j44|) are now satisfied up to the 0{e) order. Collecting terms at the O(e^) 
order, we obtain 

uS!!., + a\^{r + 2qnr^^uSl, = -2Xmill, - 2Xm!il, - {M^^Mli 
+ V"{-^{t + 2qn))wil^ + V"{ip{T + 2qn))W^Z2 



V"'{-^{t + 2qn)){wi^\T + 2qn) - u^^\t + 2qn))uf^^ 
V"'{ip{T + 2gn))(u?V + 2gn) - w^^\t + 2qn - 2q))UZ'-i, 



2n-l (55) 



r(i) w(0)^ , ^r'ff , o.^^^^r^(l) w^), 



+ y"((^(r + 2qn + 2g))([/£Vi " ^i?) + ^"("'^(^ + 2'Zn))(^£li " ^2 



where corrections ul^^ and w^^ are defined by expansion (j26p . 



To solve the linear inhomogeneous equations ()55p , the source terms have to satisfy the Fredholm 
conditions because both operators L and Lq defined by ()29p and ()30p have one-dimensional kernels. 
Therefore, we require the first equation of system (I55p to be orthogonal to <f and the second 
equation of system ([55]) to be orthogonal to 1 on [— tt, vr]. Substituting ([16]) and ([M]) to the 
orthogonality conditions and taking into account the symmetry between couplings of lattice sites 
on Z, we obtain the difference equations for {c2n-i,a2n}nez- 



Kh?C2n-l = Mi{c2n+1 + C2n-3 " 2C2n-l) + LiA{a2n - a2n-2), ^^q^ 



h?a2n = M2{a2n+2 + a2n-2 " 2a2„) + L2A(c2n+l - C2n-l), 
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where A = A^^^ and {K, Mi, M2, Li, L2) are numerical coefficients to be computed from the 
projections. In particular, we obtain 



K = {2v{t) + (fir)) (p{T)dT, 

J — 7r 

Ml = r V"{-ip{T))ip{T)z+{T + 2q)dT = r V"{ip{T))ip{T)z^{T - 2q)dT, 

J —TV J —TV 

M2 = ^ V"{^{T + 2q))y^{T + 2q)dT = ^f V" {-^{T))y+{T)dT, 
Li = -2 I y^{T)^{T)dT = 2 j y+iT)ipiT)dT, 

J —TV J —TV 

L2 = V"{^{t + 2q))v{T + 2q)dT = -i- y"(-(^(r))i;(T)dT. 

Note that the coefficients Mi and M2 need not to be computed at the diagonal terms C2n-i 
and 02n thanks to the fact that the difference equations (j56p with A = must have eigenvectors 
with equal values of {c2„,_i}nGZ and {a2n}n€Z, which correspond to the two symmetries of the 
linearized dimer system (|38p related to the symmetries ([7]) and dS]). This fact shows that the 
problem of limited smoothness of V" , which is C but not near zero, is not a serious obstacle 
in the derivation of the reduced system ([56|) . 

Difference equations (j56p give a necessary and sufficient condition to solve the linear inhomoge- 
neous equations (I55p at the ©(e^) order and to continue the perturbation expansions beyond this 
order. Before justifying this formal perturbation theory, we shall explicitly compute the coefficients 
[K, Ml, M2, Li, L2) of the difference equations ([56]) . 

Note that the system of difference equations ([56|) presents a quadratic eigenvalue problem 
with respect to the spectral parameter A. Such quadratic eigenvalue problem appear often in the 
context of spectral stability of nonlinear waves [5l [16] . 

4.4 Explicit computations of the coefficients 

We shall prove the following technical result. 

Lemma 1. Coefficients K, M2, Li, and L2 are independent of q and are given by 

K M 2 2{2,r-T'{E,){m?) 



Consequently, K > 0, whereas M2, Li, L2 < 0. On the other hand, coefficient Mi depends on q 
and is given by 

Mi = --i^{0)f+I{q), 



IT 

where 



I{q) = I{tt -q):=- (p{t)(p{t + 2q)dT, q G 

J-rT-2q 



TT 

"•2 
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To prove Lemma [H we first uniquely solve the linear inhomogeneous equations ([5T]) . ([52|) . and 
(f53]) . For equation (i5T]l . we note that a general solution is 

v{t) = -tlp{t) + 6i0(t) + b2dE^PEo{'r), r e [-7r,7r], 

where (61,62) are arbitrary constants and Be'^Eo is the derivative of the r(i?)-periodic solution 
of the nonlinear oscillator equation (jl6p with the first integral (jl7p satisfying initial conditions 
'/'b(O) = and ^e{^) = V^E' at the value of energy E = Eq, for which T{Eq) = 27r. We note the 
equation 

dEipEoi^^) = TlT'{Eo)^{±7r), (57) 

that follows from the differentiation of equation ipE{iT{E)/2) = with respect to E at E = Eq. 

To define v uniquely, we require that {ip,v)i2 = 0. Because 99 is even in r, whereas np and 
Qe^Eo are odd, we hence have 61 = and v{0) = 0. Hence v is odd in r and, in order to satisfy the 
27r-periodicity, we shall only require v{7r) = 0, which uniquely specifies the value of 62 by virtue 
of dSTD, 

dE^EoiT^) ^~T'{Eo)' 

As a result, we obtain 

v[r) = -Tipir) - j^^—:^dEipEo{T), TG[-7r,7r]. (58) 

For equation (f52|) . we can use that (/^(r) > for r G [0, vr] and (^(r) < for r G [— 7r,0]. We 
can also use the symmetry (^(tt) = —99(0). Integrating equations for y^. separately, we obtain 



y+i-T) 
y-{T) 



l+a+ip + b+dE^PEo, rG[-7r,0], 

C+if + d+dE^PEo, 'T ^ [0' 

a^if + b^dE^PEo, rG[-7r,0], 

1 + c-(p + d_dE<fEo, t"G[0, vr]. 



Continuity of y± and y± across r = defines uniquely d± = b± and c± = a± i: With this 
definition, y±{—iT) = y±(7r), whereas condition y±(— vr) = y±(vr) sets up uniquely 

2 



T'{Eo)m' 



whereas constants a± are not specified. 

To define y± uniquely, we again require that {^,y±)L'^^^ = 0. This yields the constraint on a±, 

1 2{ip,dEVEo)Ll^^ 

a± = =F: 



2m T'{Eo)m{^.^)LlJ 
As a result, we obtain 

f 1, r G [-7r,0], 
2/+(t) = a+(^(r) + b+dE^Eoir) + < ^(r) ^ ^ rn ^1 
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and 



y_(r) = a_(^(T) + b-dE^PEoi^) + 



0, Te[-7r,0], 



(60) 



where (a±,6±) are uniquely defined above. 

For equation ([53]) . we integrate separately on [— vr,0] and [0,7r] to obtain 



c+-\ip{T)r, TG[-7r,0], 
c+, T G [0,7r], 



i_(r) 



c_, r G [-7r,0], 

c„ + |(/.(r)r, rG[0,7r], 



where (c+,c_) are constants of integration and continuity of z± across r = have been used. To 
define z± uniquely, we require that {1,z±)e2 = 0. Integrating the equations above under this 
condition, we obtain: 



c+T + d+ - ip{t), TG[-7r,0], 



c_i.r — dj 



T G [0, vr] 



z.{r) 



c^T + d-, T G [— vr, 0], 

C-T — d- — ip{T), TG[0,7r], 



where {d+,d-) are constants of integration. Continuity of z± across r = uniquely sets coefficients 
d± = ibi(/3(0). Periodicity of z±{—ir) = z±{ir) is satisfied. Periodicity of z±{—ir) = z±{ir) uniquely 
defines coefficients c± = ±-(^(0). As a result, we obtain 



1 / 0(O)(2r + 7r) - 27r(/3(r), TG[-7r,0], 

r G [0,7r], 



271 \ 9i(0)(2r-7r), 



and 



^ ' 27r 



(61) 
(62) 



-V3(0)(2r + ^), TG[-^,0], 
-V3(0)(2r - vr) - 27rvi(T), TG[0,7r]. 

We can now compute the coefficients {K, Mi, M2, Li, L2) of the difference equations (|56|) . For 
coefficients K, we integrate by parts, use equations (fT6|) . pT|) . ([58]) . and obtain 



K 



ip{ip + 2v)dT 
2tt 



{if^ — 2v(p)dT 



T'{Eo) 



+ 



2tt 



dp 



T=~7V 

1 -2 1 

-if H 

2^ l + a 



T'{Eo) 

l+a 



dr 



47r2 



T'{Eo 



Because T'{Eq) < 0, we find that A' > 0. 

For Ml, we use equation (f53]l . solution ([62]) . and obtain 



Ml 



y"(-(^(T))vi(T)z+(r + 2g)(ir = / z_(r)z+(r + 2g)dT 

-TT ./ — TT 



i_(r)i+(r + 2q)dT 



(^(r)i+(r + 2g)(ir, 
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hence, the sign of Mi depends on q. Using solution ([62]) . for q G [O, ^] , we obtain 

1 



Ml 



vr 



9^(0) / Cp{T)dT 



(f{T)(p{T + 2q)dT 



TT 



:(0(O))2+/(g), liq) 



(P{t)(p{t + 2q)dT. 



TT-2q 



On the other hand, for g E vr] , we obtain 

2 



Ml 



TT 



;(9i(0))^ + /(g), /(g) := 



27r-2g 



SO that 



/(vr 



(/3(r)(/3(r — 2q)dT 



2q 



(^(T)(/3(r + 2g)(ir, 



¥3(T)¥3(T + 2g)dT = /(g), 



because the mean value of a periodic function does not depend on the limits of integration. 
For M2, we use equation ([52]) and obtain 

1 



Mo 



27r 

1 

~2^ 



y"(-(^)y+(ir 



a 
2^ 







1 • 2 

vr 7rT'(Eo)(v5(0))^ 



hence, M2 < 0. 

For Li, we use equations (fTUj) . (HZD, (pUj) . and obtain 



-2 / y^vidr = -2b- / dE^Eo^dT 

— TV 

4 



T'(i?o)0(O) 
2(2^-r(i;o)(^(0))^) 

r'(£;o)0(o) 

Because (/3(0) > and T'{Eq) < 0, we find that Li < 0. 



(dE^PEo'f - dE(PEoV>) dr + '^dEipEo 



T=0. 



For L2, we use equations (jT6]) . (|T7|) . (I58]l . and obtain 



1 



27r 



a 



— / y"(-¥;)?;dT = / (/j"-^w(ir 



27r 



1 



1 



-dE'fEo 



21^^ T'{Eo 

2TT-T'{Eo){m? _ 1 
7rr'(^o)'/i(0) 27r 

hence, -L2 < 0. 

The proof of Lemma [1] is complete. 



r=0 
^1, 
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4.5 Eigenvalues of the difference equations 



Because the coefficients {K, Mi, M2, Li, L2) of the difference equations ([56|) are independent of n, 
we can solve these equations by the discrete Fourier transform. Substituting 



C2n-1 



i2en 



(63) 



(64) 



where 9 G [0, vr] is the Fourier spectral parameter, we obtain the system of linear homogeneous 
equations, 

KA'^C = 2Mi(cos(26l) - 1)C + 2iLiAsin(6l)^, 
A'^A = 2M2(cos(26') - 1)^ + 2iL2Asin(6')C. 

A nonzero solution of system (j64p exists if and only if A is a root of the characteristic polyno- 
mial, 

D{A; 9) = KA^ + 4A2(Mi + KM2 + L1L2) sm^{e) + I6M1M2 sin''(6') = 0. (65) 

Since this equation is bi-quadratic, it has two pairs of roots for each 9 G [0, vr]. For = 0, both 
pairs are zero, which recovers the characteristic exponent A = of algebraic multiplicity of (at 
least) 4 in the linear eigenvalue problem (|44p . For a fixed 9 G (0,'7r), the two pairs of roots are 
generally nonzero, say A^ and Ag. The following result specifies their location. 

Lemma 2. There exists a ^ (O, ^) such that A^ < A| < for q G [0,(70) U (vr — qo,7r] and 
Al<0<AlforqG {qo,TT-qo). 



(66) 



To classify the nonzero roots of the characteristic polynomial (|65p . we define 

r := Mi+ KM2 + L1L2, A:= 4.KM1M2. 
The two pairs of roots are determined in Table I. 



Coefficients 


Roots 


A < 


Af < < A2 


< A < F^ F > 


A^ < A^ < 


< A < F^ F < 


A^ > A^ > 


A>F^ 


Im(A^) > 0, Im(A^) < 



Table I: Squared roots of the characteristic equation (I65p . 

Using the explicit computations of the coefficients {K, Mi, M2, Li, L2), we obtain 



F 



T'{Eo 



+ I{q), A 



64 



2{mf 



Because I{q) is symmetric about g = ^, we can restrict our consideration to the values q G [O, ^] 
and use the explicit definition from Lemma [1) 



(P{t)(p{t + 2q)dT, qe 



TT-2q 



0, 



TT 



We claim that I{q) is a positive, monotonically increasing function in [O, ^] starting with /(O) = 0. 
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Because (/3(r) = — |(^(t)|"~-^(/?(t), we realize that (/3(r) < for r G [0,7r], whereas (/3(r + 2g) > 
for r e [tt — 2g,7r]. Hence, /(g) > for any 2q G [0, vr]. Moreover, / is a continuously differentiable 
function of q, because the first derivative, 

/•TT /"TT 

I'{q) = -2 / ^{T)'^{T + 2q)dT = 2 V {t)(p{t + 2q)dT 

J n—2q J ■K—2q 



= -2a / |(/,(r)r-V(r)(^(r + 2g)dT, 

JiT~2q 

is continuous for all 2q G [0,7r]. Because 93(r) and (^(r) are odd and even with respect to r = ^, 
respectively, and (^(r) > for r G [O, ^] , we have I'{q) > for any 2q G [0, vr]. Therefore, I{q) is 
monotonically increasing from /(O) =0 to 

/ TT \ 

^\2)^~Jo ^(^)^(^ + = mT))^dT>0. 
Hence, for all q G [O, , we have F > and 

^2 . w N / w N 16 327r \ 

where A = if and only if g = 0. Therefore, only the first two lines of Table I can occur. 

For q = 0, /(O) = 0, hence Mi < 0, A > and A = F^. The second hne of Table I gives 
= < 0. All characteristic exponents are purely imaginary and degenerate, thanks to the 

explicit computations: 

Al = Al = -^sm\e). (67) 

The proof of Lemma [2] is achieved if there is go G (O, ^) such that the first line of Table I yields 
Af < < A^ for g G (go, f] and the second hne of Table H yields Af < A^ < for g G (0,go). 
Because / is a monotonically increasing function of g and A > for g = 0, the existence of 
go G (O, I) follows by continuity if A < for g = |. Since A' > and M2 < 0, we need to prove 
that Ml > for g = I or equivalently. 

Because (p is a 27r-periodic function with zero mean, Poincare inequality yields 



On the other hand, using equations (fT6|) . (fT7|) . and integration by parts, we obtain 



{^{r)YdT = -- ^{T)(p{T)dT = - |(^(r)r+Mr= ' E, 
2 J-TT 2 2 (a + 3) 

where the last equality is obtained by integrating the first invariant (jl7|) on [— 7r,7r]. Therefore, 
we obtain 

vrx ^ 27r(a + 1) ^ _ 7r(a + 1) , ^ , 2 ^ 2 , ...^ 



![-]> ' '' E = -^—-^(^(0))^ > -i^m 

(a + 3) (a + 3) tt 
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Figure 1: Coefficients F (left) and A (right) versus q. 



where the last inequality is obtained for a = | based on the fact that ~ 2.74 > 1. Therefore, 
Ml > and hence, A < for q = ^. The proof of Lemma [2] is complete. 

Numerical approximations of coefficients F and A versus q is shown on Figure [H We can see 
from the figure that the sign change of A occurs at qo ~ 0.915. 

4.6 Krein signature of eigenvalues 

Because the eigenvalue problem (j64p is symmetric with respect to reflection of 9 about ^, that 
is, sin(0) = sin(7r — 6), some roots A E C of the characteristic polynomial (|65|) produce multiple 
eigenvalues A in the linear eigenvalue problem (|44p at the 0(e) order of the asymptotic expansion 
(|47|) . To control splitting and persistence of eigenvalues A G iM_|_ with respect to perturbations, 
we shall look at the Krein signature of the 2-form a defined by (j45p . The following result allows 
us to compute a asymptotically as e — t- 0. 

Lemma 3. For every q G (0,go)j the 2-form a for every eigenvector of the linear eigenvalue 
problem |.^^[ ) generated by the perturbation expansion associated with the root A G iR+ of the 
characteristic equation k65\) is nonzero. 

Using the representation (j43p for \ = iu with oo G M+, we rewrite a in the form: 

a = 2W ^ [\U2n-l\^ + \W2n?'] + i ^ \U2n~lU2n~l " U2n~lU2n-l + W^nW^n - WinWin ■ 

Now using perturbation expansion uj = e^-\- O(e^), where A = G iM+ is a root of the charac- 
teristic equation (j65]) . and the perturbation expansions ([18]) for the eigenvector, we compute 

where 

£7^1) = 2J7 [|c2„-l|2vi2(^ + 2qn) + \a2n?\ + i(c2„-lt^2n-l - C2n,-lf/2n-l)'/i(r + Iqu) 

-i(c2„,„if72?„i - C2„_it/i;,Li)v3(r + Iqn) + i^a^Jll^^ - a^nW^^})- 
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Using representation (|Mj) . this becomes 

^n^ = 2f](|c2n-lP-E'o + |a2nP) + «(c2n-ia2n " C2„-ia2n)-E- + i{c2n-ia2n-2 " C2n-ia2n-2)-E'+ , 

where £"0 and Ej- are numerical coefficients given by 

Eo = if^ + ipi) - (pv, 
E± = ipy± - (py± - z±. 

Using exphcit computations of functions v, y±, and z± in Lemma [U we obtain 

p 27r _ 2^-r(i^o)(93(0))2 

° T'(i?o)' ^ vrr'(i^o)<^(0) ' 

and hence we have 

cr^^ = 2Q ^^I'^Sn-lP + |«2n|^^ — ^-^2 (C2n-ia2n — C2n-l«2n — C2n.-ia2n.-2 + C2n-ia2n.-2)- 

Substituting the eigenvector of the reduced eigenvalue problem (j56p in the discrete Fourier 
transform form (|63p . we obtain 

= 20 (^^2 + - 4L2 sin(0)C^ 
= ^ {n^KC^ + SvrMa sin2(e)^2^ , 

where the second equation of system (j64p has been used. Using now the first equation of system 
we obtain 

("rP = r r o'^ [^^1^2^^^ + M2(i^fi2 - 4Mi sin2(0))2] . (68) 

Note that cii^^ is independent of n, hence periodic boundary conditions are used to obtain a finite 
expression for the 2-form a. 

We consider q G (0, go) and G (0, vr), so that $7^0 and C 7^ 0. Then, cj^j^^ = if and only if 

KLiL2n^ + M2{Kn'^ - 4Mi sin2(6l))2 = 0. 

Using the explicit coefficients in Lemma [H we factorize the left hand side as follows: 

KLiL2n'^ + M2iKn'^ - 4Mi sm^{e)f = (0^ + T'{Eq)MiM2 s\n^{e)) 

(T'(Eo))2 V 4, ; T'(E„) ' ' ') ' ' 

For every q S (0, go), Mi < 0, so that the second bracket is strictly positive (recall that T'{Eq) < 0). 
Now the first bracket vanishes at 

7r(v9(0))^ 
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Substituting this constraint to the characteristic equation (j65p yields after straightforward com- 
putations: 

8Afisin'^(0) / 27r 



which is nonzero for all q G (0, go) and 9 G (0, vr). Therefore, does not vanish if g G (0, qo) and 
9 £ (0,7r). By continuity of the perturbation expansions in e, a does not vanish too. The proof of 
Lemma [3] is complete. 



Remark 8. For every q £ (0,(7o)> ^,11 roots A £ iM_|_ of the characteristic equation i65\) are divided 
into two equc 
factorization 



into two equal sets, one has > and the other one has < 0. This follows from the 



As q ^ 0, I{q) —7- and perturbation theory for double roots |g7p for q = yields 



Using the factorization formula I169\). the sign of is determined by the expression 



+ T'iEo)M,M2sm\9) = ±^sm\9)^\T'{Eo)\Iiq) - ^Eo){m? ) 

which justifies the claim for small positive q. By Lemma\^ the Krein signature of does not 
vanish for all q G (0,go) o.'iT'd 9 G (0, tt), therefore the splitting of all roots A G iR+ into two equal 
sets persists for all values of q £ (0, go). 

4.7 Proof of Theorem [2] 

To conclude the proof of Theorem [2l we develop rigorous perturbation theory in the case when 
q = ^ for some positive integers m and N such that 1 < m < A^. In this case, the linear eigenvalue 
problem (j44p can be closed at 2mN second-order differential equations subject to 27nA'^-periodic 
boundary conditions (jlSp and we are looking for AmN eigenvalues A, which are characteristic 
values of a AmN x AmN Floquet matrix. 

At e = 0, we have 2mN double Jordan blocks for A = 0. The 2mN eigenvectors are given by 
(I46p . The 2mA^-periodic boundary conditions are incorporated in the discrete Fourier transform 
§B if 

9 = ^^=9k{m,N), fc = 0,l,...,miV- 1. 
mN 

Because the characteristic equation (j65p for each 9k{m,N) returns 4 roots, we count AmN roots 
of the characteristic equation (I65p . as many as there are eigenvalues A in the linear eigenvalue 
problem (j44p . As long as the roots are non-degenerate (if A 7^ F^) and different from zero (if 
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A 7^ 0), the first-order perturbation theory predicts sphtting of A = into symmetric pairs of 
non-zero eigenvalues. The zero eigenvalue of multiplicity 4 persists and corresponds to the value 
9Q{m,N) = 0. It is associated with the symmetries of the dimer equations ^ and dH]). 

The non-zero eigenvalues are located hierarchically with respect to the values of sin^(0) for 
9 = 6k[m,N) with 1 < k < niN — 1. Because sin(0) = sin(7r — 9), every non-zero eigenvalue 
corresponding to 9k{m,N) 7^ f is double. Because all eigenvalues A E iM+ have a definite Krein 

signature by Lemma[3]and the sign of in ()68p is same for both eigenvalues with sin(0) = sin(7r — 
9), the double eigenvalues A G zM are structurally stable with respect to parameter continuations 
[3] in the sense that they split along the imaginary axis beyond the leading-order perturbation 
theory. 

Remark 9. The argument based on the Krein signature does not cover the case of double real 
eigenvalues A G ]R_|_, which may split off the real axis to the complex domain. However, both real 
and complex eigenvalues contribute to the count of unstable eigenvalues with the account of their 
multiplicities. 

It remains to address the issue that the first-order perturbation theory uses computations 
of V'" , which is not a continuous function of its argument. To deal with this issue, we use a 
renormalization technique. We note that if is a solution of the differential advance-delay 

equations (|14p given by Theorem [U then 

u^:{t) = V'\ew^{T) - u^{T)){ew^{T) - u^.{t)) 

- y"K(r) - ew^{T - 2q)){u,{T) - ew^{T - 2q)), (70) 

where the right-hand side is a continuous function of r. 
Using ([70|) . we substitute 

U2n~l = C2n-lU*{T + 2qn) + U2n-1, W^2n = >V2n, 

for an arbitrary choice of {c2n-i}neZ; into the linear eigenvalue problem (j44p and obtain: 

' U2n-i + 2A^2n-i + A^Z^sn-i = V" {ew^{T + 2qn) - u*(r + 2qn)){eW2n - l^2n-i) 

- V"{u^{t + 2qn) - ew^{T + 2qn - 2q)){U2n-i - eW2n-2), 

- (2Ati*(r 2qn) ^ }?u^{t ^ 2qn))c2n~i 

- £V"{ew^,{T + 2qn) — u^,{t + 2qn))w^{T + 2qn)c2n-i 

< - eV"{u^{T + 2qn) - ew^{T + 2qn -2q))w*{T + 2qn -2q)c2n-i-, (71) 

W2n + 2\W2n + A^Wan = eV"{u,{T + 2qn + 2q) - ew,{T + 2qn)){U2n+i - eW2n) 

- eV"{ew^{T + 2qn) - u^{t + 2qn)){eW2n -l^2n~i) 

+ eV"{u^{T + 2qn + 2q) - ew^{T + 2gn))n*(r -|- 2qn + 2g)c2n-i 
-I- eV"{ew^,{T + 2qn) — u^{t + 2qn))u^,{T + 2qn)c2n-i- 

When we repeat decompositions of the first-order perturbation theory, we write 

A = eA«+e2;^(2) + o(e2), 

U2n-i = e4lli + eHn-i+o{e^), 
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for an arbitrary choice of {a2n}ngz- Substituting this decomposition to system ([7T]) . we obtain 
equations at the 0{e) and 0{e'^) orders, which do not require computations of V'". Hence, the 
system of difference equations (j56p is justified and the splitting of the eigenvalues A at the first 
order of the perturbation theory obeys roots of the characteristic equation (I65p . Persistence of 
roots beyond the o(e^) order holds by the standard perturbation theory for isolated eigenvalues of 
the Floquet matrix. The proof of Theorem [2] is complete. 



5 Numerical Results 

We obtain numerical approximations of the periodic travelling waves (1120 in the case q = jj-, where 
is an integer, when the dimer system ^ can be closed as the following system of 2N differential 
equations: 

U2n-lit) = {eW2n{t) - U2n-l{t))% - {U2n-l{t) - £W2n-2{t))%, l<n<N (72) 
W2n{t) = e{u2n-l{t) - £W2n{t))% - e{eW2n{t) - U2n+l{t))%, - - , 

subject to the periodic boundary conditions 

n„i = U2N-1, U2N+1 = Ui, Wo = W2N, W2N+2 = W2- (73) 

The periodic travelling waves (fT2]) corresponds to 27r-periodic solutions of system (f72]l satisfying 
the reduction 

U2n+l{t) = U2n-1 (^t + , W2n+2(.t) = W2n (^t + , t G M, I < U < N. (74) 

For convenience and uniqueness, we look for an odd function ui{t) = —ui{—t) with 

ni(0) = and ni(0) > 0. (75) 

By Theorem [H the travelling wave solutions satisfying (I74p and (j75p exist uniquely at least for 
small values of e. We can continue this branch of solutions with respect to parameter e in the 
interval [0, 1] starting from the limiting solutions obtained at e = 0. 



5.1 Existence of travelling periodic wave solutions 

In order to obtain 27r-periodic traveling wave solutions to the nonlinear system ()72p . we use the 
shooting method. Our shooting parameters are given by the initial conditions 

{{u2n-l{0),U2n-l{0),W2n{0),W2n{0)}l<n<N- 

Since ui(0) = 0, this gives a set of 2N — 1 shooting parameters. However, for solutions satisfying 
the travelling wave reduction ()74p . we can use symmetries of the nonlinear system of differential 
equations (I72p to reduce the number of shooting parameters to parameters. 

For two particles (A^ = 1 or g = vr), the existence and stability problems are trivial. The exact 
solution (I23p is uniquely continued for all e € [0, 1] and matches the exact solution of the granular 
chain of two identical particles at e = 1 considered in [12]. This solution is spectrally stable 
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with respect to 2-periodic perturbations for all e G [0, 1] because the characteristic value A = 
has algebraic multiplicity four, which coincides with the total number of admissible characteristic 
values A. 

For four particles (A^ = 2 or g = the nonlinear system (I72p is written explicitly as 

Ui{t) = {£Wi{t) - Ui{t))% - {Ui{t) - £W2{t))%, 

W2{t) = e[{ui{t) - ew2{t))l - {ew2{t) - u^{t)r+, , . 

U3{t) = iew2it) - U3{t))% - (usit) - ew^{t))%, ^'"^ 
Mt) = 4iu3{t) - ewi{t))l - (ew^it) - ui{t))%. 

We are looking for 27r-periodic functions satisfying the travelling wave reduction: 

U3{t) =Ui{t + TT), w^it) =W2it + Tr). (77) 

We note that the system (I76p is invariant with respect to the following transformation: 

ui{-t) = -ui(t), W2{-t) = -Wi{t), uz{-t) = -u^{t), Wi{-t) = -W2(t). (78) 

A 27r-periodic solution of this system satisfying (j78|) must also satisfy ui(7r) = U3{tt) = and 
^2(71") = —W4{Tr). Then, the constraints of the travelling wave reduction (fTTl) yields the additional 
condition W4{'it) = 102(0). 

To approximate a solution of the initial- value problem for the nonlinear system (176p satisfying 
()78p . we only need four shooting parameters (01,02503,04) in the initial condition: 

ui(0)=0, ni(0) = oi, ^2(0)= 02, -^2(0)= 03, 
7x3(0) = 0, -113(0) = 04, ^4(0) = -02, Wi{0) = 03. 

The solution of the initial-value problem corresponds to a 27r-periodic travelling wave solution 
only if the following four conditions are satisfied: 

ui{tt) = 0, W2{it) + W4(7r) = 0, 102(0) — ^4(71) = 0, ii3('/r) = 0. (79) 

These four conditions fully specify the shooting method for the four parameters (01,02,03,04). 
Additionally, the solution of the initial-value problem must satisfy two more conditions: 

W2{it) - Wiiir) = 0, W2{0) - Wiiir) = 0, (80) 

but these additional conditions are redundant for the shooting method. We have been checked 
conditions (I80p apostoreori, after the shooting method has converged to a solution. 

We are now able to run the shooting method based on conditions ()79p . The error of this 
numerical method is composed from the error of an ODE solver and the error in finding zeros for 
the functions above. We use the built-in MATLAB function ode 113 on the interval [0, vr] as an 
ODE solver and then use the transformation (|78p to extend the solutions to the interval [— vr, vr] 
or [0,27r]. 

Figure [2] (top left) shows three solution branches obtained by the shooting method by plotting 
^2(0) versus e. The first solution branch (labeled as branch 1) exists for all e G [0, 1] and is shown 
on the top right panel for e = 1. This branch coincides with the exact solution (|22p . The error 
in the supremum norm between the numerical and exact solutions ||ni — can be found in 

Table II. 
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Figure 2: Travelling wave solutions for N = 2: the solution of the dimer chain continued from 
e = to e = 1 (top right) and two solutions of the monomer chain at e = 1 (bottom left and 
right). The top left panel shows the value of W2{0) for all three solutions branches versus e. 



AbsTol of Shooting Method 


AbsTol of ODE solver 




0(10-^^) 


0(10-^^) 
0(10^10) 


4.5 X 10"^'' 
3.0 X 10"^^ 


o(io-«) 


O(lO^i^) 

o(io-io) 


4.5 X lO"^"* 
3.0 X 10-11 



Table II: Error between numerical and exact solutions for branch 1. 

We can see from the top left panel of Figure [2] that a pitchfork bifurcation occurs at e = eo ~ 
0.72 and results in the appearance of two symmetrically reflected branches (labeled as branches 2 
and 2'). These branches with 1^2(0) 7^ extend to e = 1 (bottom panels) to recover two travelling 
wave solutions of the monomer chain ([6|). The solution of branch 2 satisfies the travelling wave 
reduction Un+i{t) = Un [t + ^) and was previously approximated numerically by James [12]. The 
other solution of branch 2' satisfies the travelling wave reduction Un+i{t) = Un{t — §) and was 
previously obtained numerically by Starosvetsky and Vakakis |22j . 

For N = 2 {q = ^), the solution of branch 2' given by {u2n-i5 w'2n}ne{i,2} is obtained from the 
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solution of branch 2 given by {u2n-i,W2n}ne{i,2}j by means of the symmetry 

Ui{t) = -U3{t), W2{t) = -W2{t), Usit) = -Ui{t), Wi{t) = -Wi{t), 



(81) 



which holds for any e > 0. (Of course, both solutions 2 and 2' exist only for e S {eq, 1] because 
of the pitchfork bifurcation at e = Eq ~ 0.72.) The solution of branch 1 is the invariant reduction 
U2n-i = U2n~-i, W2n = W2n with respect to the symmetry ([HT]) so that it satisfies W2{t) = Wiit) = 
for all t. 

For six particles {N = 3 or q = ^), the nonlinear system (j72p is written explicitly as 



f iliit) = i^Mt) - - (Mt) - ew2{t))l, 

W2{t) = e[{ui{t) - ew2{t))% - {ew2{t) - u^it))^, 
uz{t) = {ew2{t) - u^{t))% - (u3(t) - ewi{t))%, 
Mt) = ^iiMt) - emit))% - iew4{t) - U5it))%, 

[ iveit) = e[{u5{t) - ewe{t))l - (eweit) - ni(t))° . 
We are looking for 27r-periodic functions satisfying the travelling wave reduction: 



(82) 



Unit) 



.3(^+1) 



Ui[t + 



Air 



2tt\ 



W2[t + 



We note that the system ()82p is invariant with respect to the following transformation: 

Ui{-t) = -Ui{t), W2{-t) = -WQ{t), U3{-t) = -U5{t), W^{-t) = -Wi{t). 



(83) 



(84) 



A 27r-periodic solution of this system satisfying (184p must also satisfy ui{tt) = W4{tt) = 0, W2{7t) = 
—wq{tt), and U3(7r) = —U5{tt). Then, the constraints of the travelling wave reduction (j83p yield 
the conditions ^3(71) = ui (^) and W4{Tr) = W2 (f). 

To approximate a solution of the initial- value problem for the nonlinear system (j82p satisfying 
(I84p . we only need six shooting parameters (01,02,03,04,05,06) in the initial condition: 

iii(0) = 0, ni(0) = oi, ^2(0) = 02, W2{0) = as, 
^3(0)= 04, ^3(0) = 05, w;4(0) = 0, 'u;4(0) = 06, 
UbiO) = -04, -115(0) = 05, W6{0) = -02, ^6(0) = 03. 

This solution corresponds to a 27r-periodic travelling wave solution only if it satisfies the following 
six conditions: 

ui(7r)=0, W2{'rr) + WQ{7r) = 0, W3(7r) + 7/5(71) = 0, 
(f) -U3(vr) =0 W2 {^) - Wiiir) = 0, w^iir) = 0. 

The six conditions determines the shooting method for the six parameters (01,02,03,04,05,06). 
Additional conditions, 

u;2(7r) - 'u;6(7r) = 0, us^n) - U5{it) = 0, ui (^'^^ - usiir) = 0, W2 (^^^ - mi-^) = 0, 
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Figure 3: Travelling wave solutions for = 3: the solution of branch 1 is continued from e = to 
e = 1 (top right) and the solution of branch 2 is continued from e = 1 (bottom left) to e = 0.985 
(bottom right). The top left panel shows the value of W2{0) for solution branches 1 and 2 versus 
e. 

are to be checked aposteriori, after the shooting method has converged to a solution. 

Figure [3] (top left) shows two solution branches obtained by the shooting method. Again, 
w;2(0) is plotted versus e. Branch 1 is continued from e = to e = 1 (top right) without any 
pitchfork bifurcation in e € (0, 1). Branch 2 is continued from e = I (bottom left) starting with 
a numerical solution of the monomer chain ([6]) satisfying the reduction Un+i{t) = C/„ (t + ^) to 
e = 0.985 (bottom right), where the branch disappears from the radars of our shooting method. 
We have not been able so far to detect numerically any other branch of travelling wave solutions 
near branch 2 for e = 0.985, hence the nature of this bifurcation will remain opened for further 
studies. 

We use the same technique for = 4 and show similar results on Figure |4l Branch 1 is 
uniquely continued from e = to e = 1 (top right), whereas branch 2 is continued from e = 1 
(bottom left) starting with a numerical solution of the monomer chain ([6]) satisfying the reduction 
Un+i{t) = Un [t + j) to e = 0.9 (bottom right), where the branch terminates. 
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Figure 4: Travelling wave solutions for = 4: the solution of branch 1 continued from e = 
to e = 1 (top right) and the solution of branch 1 continued from e = 1 (bottom left) to e = 0.9 
(bottom right). The top left panel shows the value of W2{0) for solution branches 1 and 2 versus 
e. 

5.2 Stability of travelling periodic wave solutions 

To determine stability of the different branches of travelling periodic wave solutions of the granular 
dimer chains (jl]), we compute Floquet multipliers of the monodromy matrix for the linearized 
system (j38p . To do this, we use the travelling wave solution obtained with the shooting method 
and the MATLAB function ode 113 to compute the fundamental matrix solution of the linearized 
system ()38p on the interval [0,27r]. 

By Theorem [21 the travelling waves of branch 1 for N = 2 {q = ^) are unstable for small values 
of £. Figure [5] (top) shows real and imaginary parts of the characteristic exponents associated with 
branch 1 for all values of £ in [0, 1]. Only positive values of Re(A) and Im(A) are shown, moreover, 
Im(A) G [O, ^] because of 1-periodicity of the characteristic exponents along the imaginary axis. 

Thanks to the periodic boundary conditions, the system of linearized equations ()39p for N = 2 
is closed at 4 second-order linearized equations, which have 8 characteristic exponents as follows. 
The exponent A = has multiplicity 4 for small positive e, and two pairs of nonzero exponents 
(one is real, the other one is purely imaginary) bifurcate according to the roots of the characteristic 
equation ()65p for 9 = ^. These asymptotic approximations are shown on the top panels of Figure 
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Figure 5: Real (left) and imaginary (right) parts of the characteristic exponents A versus e for 
N = 2 for branch 1 (top) and branch 2 (bottom). 

[5] by solid curves, in excellent agreement with the numerical data. We can see that the unstable 
real A persist for all values of e in [0, 1]. The pitchfork bifurcation at e = eg ~ 0.72 in Figure [2] 
(top left) corresponds to the coalescence of the pair of purely imaginary characteristic exponents 
on Figure [5] (top right) and appearance of a new pair of real characteristic exponents for e > Eq 
on Figure [5] (top left). Therefore, the branch continued from e = is unstable for all e G [0, 1]. 

Bottom panels on Figure [5] shows real and imaginary parts of the characteristic exponents 
associated with branch 2 (same for 2' by symmetry) for all values of e in [sq, 1]. We can see that 
these travelling waves are spectrally stable near e = 1 in agreement with the numerical results of 
James jl2j . When e is decreased, these travelling waves lose spectral stability near e = ei w 0.86 
because of coallescence of the pair of purely imaginary characteristic exponents and appearance 
of a new pair of real characteristic exponents for e < £i. The two solution branches disappear as 
a result of the pitchfork bifurcation at e = eg ~ 0.72, which is again induced by the coalescence of 
the second pair of purely imaginary characteristic exponents. 

For N = 3 {q = | ) , the system of linearized equations (p9]) is closed at 6 second-order linearized 
equations. Besides the characteristic exponent A = of multiplicity four, we have 8 nonzero 
characteristic exponents A. The characteristic equation (j65|) with = ^ and 9 = ^ predicts a 
double pair of real A and a double pair of purely imaginary A. Figure [6] (top) shows Re(A) (left) 
and Im(A) (right) for solutions of branch 1. The double pair of purely imaginary A split along 
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Figure 6: Real (left) and imaginary (right) parts of the characteristic exponents A versus e for 
= 3 for branch 1 (top) and branch 2 (bottom). 

the imaginary axis for small e > 0. On the other hand, the double pair of real A splits along the 
transverse direction and results in occurrence of a quartet of complex- valued A for small e > 0. 
These complex characteristic exponents approach the imaginary axis at e = ei ~ 0.43 (Neimark- 
Sacker bifurcation) and then split along the imaginary axis as two pairs of purely imaginary A for 
£ > £\. At the same time, one pair of of the purely imaginary A continued from e = approaches 
the line ±| (corresponding to the Floquet multiplier at — 1) at e = £2 ~ 0.72 (period-doubling 
bifurcation) and splits along the negative real axis. In summary, the periodic travelling wave of 
branch 1 for = 3 is stable for e € (ei, £2) but unstable near e = and e = 1. 

Figure [6] (bottom) shows Re(A) (left) and Im(A) (right) for solutions of branch 2 that exists 
only for e S [e*, 1], where e* w 0.985. All four pairs of the characteristic exponents A are purely 
imaginary near e = 1 that corresponds to the numerical results for stability of travelling waves in 
monomer chains in [12]. Two pairs coalesce at e ~ 0.995 resulting in the complex characteristic 
exponents (Neimark-Sacker bifurcation). One more pair crosses the line ib| for e ~ 0.989 resulting 
in the negative characteristic exponents (period-doubling bifurcation). The last remaining pair of 
purely imaginary A crosses zero near e = ~ 0.985 that indicates that termination of branch 2 is 
related to a local bifurcation. However, we are not able to identify numerically any other branch 
of travelling wave solutions in the neighborhood of branch 2 for e ~ e^,. 

Recall that the coefficient Mi changes sign at g ~ 0.915, as seen in Figure [TJ Therefore, for 
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Figure 7: Real (left) and imaginary (right) parts of the characteristic exponents A versus e for 
= 4 for branch 1 (top) and branch 2 (bottom). 

> 4, the characteristic equation (j65p for any values of 9 predicts pairs of purely imaginary A 
only. This is illustrated on the top panel of Figure [7] for N = A [q = We can see that all 
double pairs of purely imaginary A split along the imaginary axis for small e > and that the 
periodic travelling waves of branch 1 remain stable for all e G [0, 1]. The figure also illustrate the 
validity of asymptotic approximations obtained from roots of the characteristic equation (|65|) . 

It is interesting that Figure [7] shows safe coalescence of characteristic exponents for larger 
values of e. Recall from Remark [S] that the characteristic exponents have opposite Krein signature 
for small values of e in such a way that larger exponents on Figure [7] have negative Krein signature 
fj and smaller exponents have positive Krein signature a. It is typical to observe instabilities after 
coalescence of two purely imaginary eigenvalues of the opposite Krein signature [T7] but this only 
happens when the double eigenvalue at the coalescence point is not semi-simple. When the double 
eigenvalue is semi-simple, the coalescence does not introduce any instabilities [3]. This is precisely 
what we observe on Figure [71 After coalescence for larger values of e, the purely imaginary 
characteristic exponents A reappear as simple exponents of the opposite Krein signature and the 
exponents with positive Krein signatures are now above the ones with negative Krein signatures. 

Figure [7] (bottom) shows Re(A) (left) and Im(A) (right) for solutions of branch 2 that exists 
only for e G [£*)!]) where e* ~ 0.90. Besides the pairs of purely imaginary characteristic exponents 
A, there exists one pair of real exponents A near e = 1 that corresponds to the numerical results for 
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Figure 8: Imaginary parts of the characteristic exponents A versus e for N = 5 (left) and = 6 
(right). The real part of all the exponents is zero. 

instability of travelling waves in monomer chains in [12]. For smaller values of e, more instabilities 
arise for the solutions of branch 2 because of various bifurcations of pairs of purely imaginary 
exponents A. 

Finally, Figure [8] illustrate the stability of solutions of branch 1 for N = 5 (left) and = 6 
(right). Not only the double pairs of purely imaginary A split safely along the imaginary axis 
for small e > 0, various coalescence of purely imaginary exponents A of opposite Krein signature 
never result in occurrence of complex exponents A. The solutions of branch 1 remain stable for all 
ee[0, 1]. 

5.3 Stability of the uniform periodic oscillations 

The periodic solution with q = (which is no longer a traveling wave but a uniform oscillation 
of all sites of the dimer) is given by the exact solution (j23p . Spectral stability of this solution is 
obtained from the system of linearized equations ()40p . Using the boundary conditions 

U2n+1 = e^'^U2n-l, W2n+2 = e^'^ W2n, n G Z, 

where 6 € [0, vr] is a continuous parameter, we obtain the system of two closed second-order 
equations, 

w + ^\v>\''-'w = j^{V"i-ip) + V"i^)e^^^')u. ^ > 

The characteristic equation (|65|) for q = Q predicts a double pair ()67p of purely imaginary A for 
any 6 G (0, vr). We confirm here numerically that the double pair is preserved for all E G [0, 1]. 

Figure [9] shows the imaginary part of the characteristic exponents A of the linearized system 
(jSSp for ^ = f (left) and 6 = ^ (right). Similar results are obtained for other values of 0. Therefore, 
the periodic solution with q = Q remains stable for all values of e G [0, 1]. 

The pattern on Figure [9] suggests a hidden symmetry in this case. Suppose \e is a characteristic 
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Figure 9: Imaginary parts of the characteristic exponents A versus e for = ^ (left) and = j 
(right). The real part of all the exponents is zero. 



exponent of the system (j85p for the eigenvector 



u 




' Ue{t) ' 


w 




_ We{t) 



(86) 



where UQ{t) and We{t) are 27r-periodic and the subscript 6 indicates that the system (jSSh depends 
explicitly on 0. Recall that the unperturbed solution satisfies the symmetry ip{t + vr) = —f{t) for 



all t. Using this symmetry and the trivial identity e 



1, we can verify that there is another 



solution of the system ()85p with the same 9 for the characteristic exponent A,, 



u 




w 





32^^VF^_e(t + 7r) 



(87) 



Prom the symmetry of roots ()67|) and the corresponding characteristic exponents, we have Xq = 
X-jT-e- The eigenvectors (j86p and (j87p are generally linearly independent and coexist for the same 
value of A = Ag = A,r-e- This argument explains the double degeneracy of characteristic exponents 
A for the case q = for all values of e G [0, 1]. 



6 Conclusion 

We have studied periodic travelling waves in granular dimer chains by continuing these solutions 
from the anti-continuum limit, when the mass ratio between the light and heavy beads is zero. We 
have shown that every limiting periodic wave is uniquely continued for small mass ratio parameters. 
Although the vector fields of the granular dimer chain equations are not smooth, we can still use the 
implicit function theorem to guarantee that the continuation is with respect to the mass ratio 
parameter. We have also used rigorous perturbation theory to compute characteristic exponents 
in the linearized stability problem. From this theory, we have seen that the periodic waves with 
the wavelength larger than a certain critical value are spectrally stable for small mass ratios. 

Numerical computations are developed to show that the stability of these periodic waves with 
larger wavelengths extends all way to the limit of equal mass ratio. On the other hand, we have 
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also computed periodic travelling waves that are continued from solutions of the granular monomer 
chains at the equal mass ratio, their spectral stability, and their terminations for smaller mass 
ratios. 

Among open problems, we have not clarified the nature of bifurcation, where the solutions of 
branch 2 terminate at a G (0, 1) for = 3, 4. We have not been able to find another solution 
nearby for e ^ e^,. Safe coalescence of purely imaginary characteristic exponents A of opposite 
Krein signatures is also amazing and we have not been able to explain the hidden symmetry that 
would explain why the eigenvalues at the coalescence point remain semi-simple. These problems 
as well as analysis of the periodic travelling wave solutions for other values of q will wait for further 
studies. 
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